Decay times in turnover statistics of single enzymes 
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The first passage times for enzymatic turnovers in non-equilibrium steady state display a statistical 
symmetry property related to non-equilibrium fluctuation theorems, that makes it possible to extract 
the chemical driving force from single molecule trajectories in non-equilibrium steady state. Below, 
we show that the number of decay constants needed to describe the first passage time distribution 
of this system is not equal to the number of states in the first passage problem, as one would 
generally expect. Instead, the structure of the kinetic mechanism makes half of the decay times 
vanish identically from the turnover time distribution. The terms that cancel out correspond to the 
eigenvalues of a certain sub-matrix of the master equation matrix for the first exit time problem. 
We discuss how these results make modeling and data analysis easier for such systems, and how the 
turnovers can be measured. 



PACS numbers: 87.10.-c,82.39.-k,02.50.-r,05.40.-a 



Enzymes are vital to most biochemical reactions, to 
increase reaction speed and as active components in cel- 
lular regulatory networks. Observations of the fluctua- 
tions on the single molecule level can lead to new in- 
sights into enzymatic mechanisms, by revealing more de- 
tailed information than ensemble averages measured in 
bulk experiments 0, [2, 0, 0, IE IE 0)- The rapid devel- 
opment in single molecule techniques has made it possi- 
ble to directly observe turnover events of single enzymes 
in many systems (see e.g., Ref. @ and refs. therein). 
This motivates continuing theoretical interest in stochas- 
tic kinetics. For example, it was recently shown how 
non-equilibrium fluctuation theorems make it possible to 
extract the chemical driving force from turnover traces 
of single enzymes 0, 

Here we consider the statistical properties of reversible 
enzymatic turnovers, and derive another useful property 
of the turnover times. Their distribution is a sum of 
exponentially decaying terms, and the number of terms 
is usually expected to reflect the number of states in 
the underlying first passage problem. We show that the 
number of terms in the actual distributions are only half 
of the expected number, due to the periodicity of the 
problem. Moreover, we discuss an earlier suggestion [1] 
of how to detect turnover events, and conclude that it 
does not correspond to the first passage problem for 
turnover times. We have previously addressed an analo- 
gous issue for stepping motor proteins, and shown that 
it can lead to systematic misinterpretations of experi- 
mental data[E0|- O ur results point to modifications in 
previously suggested experiments, and also simplify the- 
oretical analysis of turnover time distributions. 

In the next section, we introduce our model and the re- 
sults. After that, we discuss how turnover times can be 
detected in reversible single molecule experiments. We 
then derive our main result, and finally discuss some im- 



plications. 

Model. Following Qian and Xie 0], we start with 
a simple sequential kinetic model of an enzyme reac- 
tion, sketched in Fig. [lja), where a substrate A is 
converted to a product B through several intermediate 
states. The overall concentrations of substrate and prod- 
uct molecules are assumed to be kept constant, so that a 
non-equilibrium steady state is maintained. 

As sketched in Fig. [Ijb), a ± turnover is defined as 
the first arrival in one of the empty states E±„, after 
start in state Eo. If a turnover is completed at t = 0, 
the integrated turnover times w±(t) is the probability 
that the next turnover is a ±, and occurs at time t or 
earlier. Thus, the turnovers are equivalent to the cy- 
cle completion events associated with the work of Hill 
Microscopic reversibility leads to a symmetry prop- 
erty for the forward (+) and backward (— ) turnover 
timesQH, namely, w+(t) = e Atl/kBT w-(t) @, Here, 



Afi 



hi 



miom2i .. .m n 
rn imi2...m„ 



and rriij is the rate of the transition Ej 



is the chemical driving force, 
Ei. 



The turnover time distributions arc of the general form 



w±(t) 



N 
k=l 



(1) 



The characteristic decay times Tfc = — 1/Afc and prefac- 
tors a k depend on the transition rates and topology of 
the underlying kinetic mechanism. Hence, this mech- 
anism can be studied by fitting theoretically predicted 
distributions to experimental data. 

The underlying first passage problem is governed by a 
system of linear master equations [TT1 | , one equation for 
each state from which the systems escapes. Since one 
generally expects a matrix of dimension ./V to have N 
eigenvalues, a simple and common way to estimate the 
number of states is to count how many exponential terms 
are needed to fit the first passage time distribution. As 
illustrated in Fig. [T] the turnover events correspond to 



escape events from states Ei 



Eo 



■ ■ ,Et; 



Hence, 
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the number of states in the first passage problem is N 
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Figure 1: (a) A multi-step enzyme reaction converting A to B. 
The turnover time distributions W±(t) are given by the first 
exit time problem in (b). The enzyme starts in state Eo and is 
absorbed in E„ or E_„. The rate constants mio and m n _i,„ 
are pseudo-first-order, i.e., proportional to the concentrations 
of A and B respectively. In the more general scheme studied 
in Ref. 0], arbitrary transitions within a cycle are allowed. 
An example of this, with n = 5 different enzymatic states, 
is illustrated in (c). By periodicity, state Efe is equivalent to 
Efe + „, and rriij = mi+ n ,j+n- 



2n — 1 in this case, where n is the number of intermediate 
states of the enzyme-substrate complex. As shown below, 
the structure of this first passage problem, as well as the 
more general one studied by Wang and Qian Q, makes 
n — 1 of the coefficients in w± (t) vanish. This leaves 
only 7i terms in the distribution, i.e., the above estimate 
fails by a factor two. Before we derive this, we discuss 
how turnover times can be detected in reversible single 
molecule experiments. 

Detecting enzymatic turnovers. Single enzyme exper- 
iments using fluorescence techniques often probe the 
state of a enzyme-substrate complexpjl, [T^, [lj, Il5j . but 
do not report directly on the number of turnovers. In 
our example model, a realistic possibility is that the 
empty states (Eo, E-t„, . . . ) can be experimentally dis- 
tinguished from the other states, but not from each other. 
If the product concentration is kept very low, it is safe to 
assume that each departure from an empty state starts a 
new forward turnover. However, detecting individual for- 
ward and backward turnovers in conditions where back- 
ward turnovers are possible is more complicated, as the 
following discussion will show. 

In their proposal to measure A/i directly from turnover 
traces, Qian and Xie [J] suggested that individual 
turnover times could be measured by monitoring the net 
number vb[€) of product molecules, as they are released 
and absorbed by the reaction E„_i ^ E„. However, this 
is equivalent to monitoring the position of a processive 
motor protein. A closer examination reveals that this sit- 
uation corresponds to a different first passage problem, 
with quite different statistical properties. This discrep- 
ancy can lead to large systematic errors in the estimate 



for A/i[|,[3. 

To see why that measurement will not detect 
turnovers, note that the turnover event starts in state 
Eo, and finishes when an enzymatic cycle is completed, 
i.e. .when either E„ or E_„ is reached for the first time 
@i M 3- However, vb(£) docs not change during the 
reaction Ei_„ — > E_„. Therefore, backward turnovers 
cannot be detected by only monitoring changes in vq. 

The attractive statistical properties of turnover 
timesgl, motivate a consideration of how they could 
be measured, using slightly different experimental se- 
tups. One possibility would be to monitor both sub- 
strate and product molecules, but one could also imag- 
ine various setups involving fluorescence techniques with 
multiple fluorescence levels, in the spirit of the experi- 
ment that demonstrated bi-directional rotation in ATP 
synthase [l6j]. To summarize, it is important to make sure 
that the theoretical first passage problem describes the 
actual experimental situation. 

Number of decay times in w± (t) . We now derive our 
main result, i.e., that the number of exponential terms 
in the turnover time distributions w± (t) are not given by 
the number of states in the first passage problem, In — 1, 
as one might expect Q. Instead, w±(t) only contains n 
terms. The decay constants that drop out are the eigen- 
values of a certain sub-matrix of the master equation 
matrix for the first exit time problem. This result may 
simplify practical calculations considerably. 

The turnover times studied by Qian and Xie [3] are the 
solutions of the first exit problem illustrated in Fig. []Jb) . 
If we label the states —n, — n+ 1, . . . , n — 1, n, the system 
starts in state and is absorbed in states ±n. Let qk(t) be 
the probability of being in state k a time t after starting 
in state 0. The qk{t) are governed by the master equation 

d t qi[t) = ^2 (///,, 1 ~ mji<li{t)\ , —n <i < n, (2) 

with initial condition qj(0) — Sjfi. Since in are ab- 
sorbing states, the integrated turnover time distribution 
functions are given by 

w±(t) = q± n (t) = J dq± ™ t 1 dt= ^2 m± n .j J qj(t)dt. 

j—l—n ^ 

(3) 

Note that w±(t) are not normalized to unity. Instead, 
the fractions p± of ± turnovers are given by p± = 
lim^oo w± (t) . Introducing the matrix M and vector q(t) 
with elements 

n-l 

Mij = m i: j - Sij ^2 m iki ~n <i,j < n i (4) 

k— 1 — n 

fit) = [q- n+1 (t),q- n+2 {t), g„_i(i)] T , (5) 

the solution of Eq. ([2]) can be written q(t) = e tM g(0). 

We restrict our attention to the generic situation where 
M can be diagonalized [l7| . In this case, q(t) can be ex- 
pressed in terms of right and left eigenvectors aS ' and 



3 



&( fc ) of M, namely 



2m- 1 



(0) 



fc=i 



Note that the eigenvalues Afc need not all be distinct [18j| . 
Among the 2n— 1 terms, we look for left eigenvectors 
that are orthogonal to the initial condition <f(0). Those 
terms drop out of Eq. ([6]), and hence from w±(t) as well. 
We introduce 



JJ+ = R-n,o, 
V + = [mo,l-n, 
tf*0 _ r;,( fe ) 



,m-ifll ,y- 
,mo,-iJ ,w- 



["11,0, ■ 

Kii ■ 

(k) 



,m n -i,o\ 



, mo«-i 



h Cfc) 



T 

(7) 



and take Y as the (n — 1) x (n 



1) matrix with el- 
ements Yy = Mjj for < i,j < n, i.e., the master 
equation matrix for the first exit problem from states 
1, 2, ... ,n — 1. Using the periodicity of the transition 
rates, m^j = mj+ n j+ n , and a 'bottleneck' property of 



state 0, m 



i,j-n 



in i 



for < i,j < n, the left 



eigenvalue problem for M can be written 



Y 

^ M 00 4 
y+ Y T 





~b {k) ' 




b (k) 







= \ k b (k \ (8) 



where is the (n — 1) x (n — 1) zero matrix. This struc- 
ture of M also holds for the more general turnover time 
problem studied in Ref. [|| and illustrated in Fig. QJc). 
In the special case studied in Ref. 1 (Fig. Bib)), M and 
Y are tridiagonal, in which case y± , v± have only one 

non-zero element each. Setting &w ■ <f(0) 
Eq. iJHJ) gives 



Y T b {k) 



Xkb 




tfe) 



Y T bf = A fc fe 



+ ■ 



in 

(9) 
(10) 
(11) 



The solutions are given by the eigenvalues of Y T , which 
are also eigenvalues of M. Since these have equal al- 
gebraic and geometric multiplicity by assumption, there 
are n — 1 solutions, corresponding to terms that do not 
contribute to the turnover time distributions in Eq. (J6]). 
Hence, w±(t) contains at most n exponential terms: 
those where Afc is an eigenvalue of M, but not of Y. 

Turnover time distribution for sequential models. An- 
alytical expressions for the turnover time distributions 
is useful for efficient parameter extraction. Our result 
for the number of exponential terms in w± (t) makes the 
derivation of such expressions easier, and extends the 
range of system sizes that can be treated analytically. 

To illustrate this, we analyze the sequential model in 
Fig. [Ua), using the ansatz 



w±(t) =p±{l + a 1 e Mt + ...+ 



a n e 



), (12) 



with cik = (X^ /p± , together with the initial conditions for 
w±{t). As shown by Qian and Xie Q (their Ref. [24]), 
the sequential models satisfy p± w±(Q) = 1 + J2k a k = 0> 
and pg 1 d^"w±{0) = J2k x T a k = for 1 < m < n - 1. 
This leads to a Vandermonde type system of equations, 



(13) 



Solving with Cramer's rule for the normalized distribu- 
tion p± 1 w± (t), we get 



1 1 . 


. 1 




Oil 




" -1 " 


Ai A2 


A n 




a 2 







\tl— 1 \7l— 1 
A l A 2 


\7l— 1 












P±w±(t) =! + (-!)" 



n 



A, 



At, — A r 



(14) 



This is the distribution of a sum N of indepen- 
dent exponential random variables with mean values 
I Ai j , I A2 1 1 , • ■ ■ , |A„| _1 . The reduced number of un- 
known coefficients olu simplifies the analytical compu- 
tation significantly, especially so for sequential models, 
where J^ fc A™afc 7^ for m > n 

Analytical calculation of the eigenvalues Afc means 
finding the roots of a characteristic polynomial. (For 
large systems, the eigenvalues must be found numeri- 
cally In these cases, root-finding in characteristic poly- 
nomial is usually not the best method[lj|.) Since the 
n — 1 non-contributing time constants are the eigenval- 
ues of Y, these can be removed from the eigenvalue equa- 
tion in advance, hence reducing the problem from root- 
finding in the characteristic polynomial of M, Pm(A) = 
dct(M — AI), which has degree 2n — 1, to root-finding in 
the polynomial Pm(A)/Py(A), which has degree n. This 
makes it feasible to compute decay constants analytically 
for larger systems. 

Conclusion We have demonstrated that enzymatic 
turnover times constitute a counterexample to the expec- 
tation that the number of states in a first passage time 
problem is equal to the number of exponential terms in 
the first passage time distribution. Instead, the number 
of terms is in this case equal to the number of states per 
cycle. This number is an important characteristic of the 
kinetic mechanism of an enzyme, and our results make 
it possible to estimate it correctly from time series of 
turnover times. 

Furthermore, our results make it easier to derive 
(semi)analytical expressions for the turnover time dis- 
tributions, thus simplifying modeling and data analysis. 
The approach demonstrated above for a sequential model 
also works for large systems, if the eigenvalues are com- 
puted numerically. 

Finally, we have supplied an important clarification 
to an earlier suggestion Q on how to detect turnover 
events. This should make our results, together with ear- 
lier predictions [1, 0, [1], into useful analysis tools for fu- 
ture experiments, and a further reason to study enzy- 
matic turnovers under reversible conditions. 
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